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The basic box and phase-space integrals needed to compute at second order the three-jet decay 
rate of the Z-boson into massive quarks are presented in this paper. Dimensional regularization is 
used to regularize the infrared divergences that appear in intermediate steps. Finally, the cancella- 
tion of these divergences among the virtual and the real contributions is showed explicitly. 
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Oh 1 I. INTRODUCTION 

<: 

■^j- \ In higher order QCD calculations at LEP the effects of quark masses are always screened by the ratio m 2 q / 's, where 
m q is the mass of the quark and s is the center of mass energy. At LEP, even for the heaviest produced quark, 
the bottom quark, we have (mb/mz) 2 _~ 0.003 and, therefore, for inclusive observables, such as total hadronic cross 
section, or the decay width of Z — > bb, the mass effects are almost negligible. The situation is different for more 
exclusive observables. It has been stressed ]lj that quark mass effects could be enhanced in jet production cross 
sections. In this case observables depend also on a new variable, y c , the jet-resolution parameter that defines the 
£f*^ • multiplicity of the jets. In fact, these effects have already been seen However, it has also been shown in 0] 
that second-order QCD predictions are necessary in order to perform a meaningful comparison of the theory with 
■ experimental results and to possibly extract the value of the quark mass from such measurements. 
0\ ' In previous papers H we have presented the results of the calculation of some three-jet fractions at LEP including 
next-to-leading (NLO) corrections and keeping the full dependence on the 6-quark mass. In this paper we focus on 
O ' several technical details that appear in the calculation. Our results for some IR divergent integrals with massive 
■ , quarks are new and could be used in many other calculations of this type. 

The main difficulty of the NLO calculation is the appearance, in addition to ultraviolet (UV) divergences, of infrared 
(IR) singularities because gluons are massless. The Bloch-Nordsieck (BN) and Kinoshita-Lee-Nauenberg (KLN) 
theorems Q assure that jet cross-sections are infrared finite and free from collinear divergences. The IR singularities 
of the NLO one-loop Feynman diagrams cancel against the IR divergences that appear when the differential cross- 
section for four-parton production is integrated over the region of phase-space where either one gluon is soft or 
two gluons are collinear. Nevertheless, infrared and collinear divergences appear in intermediate steps and should 
be treated properly. We used dimensional regularization to regularize both UV and IR divergences || because it 
preserves the QCD Ward identities. 

For massive quarks, the gluon-quark collinear singularities, that appear in the massless case, are softened into 
logarithms of the quark mass. The number of IR singularities is, therefore, smaller. However, quark masses complicate 
extremely the loop and phase-space integrals. In particular, diagrams containing gluon-gluon collinear divergences are 
much harder to calculate than in the massless case. Furthermore, infrared singularities appear always at the border 
of the integration region. For virtual contributions, at the border of the one unit side box defined by the Feynman 
parameters. For contributions coming from the emission of real gluons, at the border of the phase-space. Thus, any 
simplification is impossible, the full calculation has to be performed in arbitrary D = 4 — 2e dimensions and expansions 
for e — > are allowed only at the end. 

Several methods of analytical cancellation of infrared singularities have been developed in the past. The most 
popular are the phase-space slicing method |(| and the subtraction method We work in the context of the 

so-called phase-space slicing method. In this case the analytical integration over a thin slice at the border of the 
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phase-space is performed in D-dimensions and the result is added to the virtual corrections. The integration over the 
rest of the phase-space, which is finite, is done numerically for D = 4. 

The structure of this paper is as follows. In section |l| we consider the most complicated, IR divergent one- loop 
box integrals that appear in the calculation of the three-jet decay rate of the Z-boson into massive quarks at NLO. 
We start this section by presenting the calculation of some simpler three-point scalar integrals, which control the IR 
structure of the box diagrams. One of the box integrals involves the three-gluon vertex and was only known in the 
case of massless quarks p0|] . We have calculated it also for the massive case by using dimensional regularization. This 
is one of the main contributions of this paper. 

Section [in] is devoted to the calculation of the basic phase-space integrals of the Z-decay into four partons wich 
lead to IR (collinear) divergent contributions to the three-jet fractions. We used the phase-space slicing method. 
When one of the gluons is soft, the calculation is performed by imposing a theoretical cut on the energy of the gluon, 
£3 < w. The w-cvit should be small enough to allow for expansions for small w in the analytical calculations and 
should not interfere with the experimental cuts. On the other hand, it should be large enough to avoid instabilities 
in the numerical integrations. For the transition amplitudes containing gluon-gluon collinear divergences we perform 
the cut on 2/34 = 2(p3 • Pa)/s, the scalar product of the momenta of the two gluons normalized to the center of mass 
energy. These contributions have to cancel exactly the appropriate infrared poles of the virtual contributions. We 
show explicitly how these basic integrals share the same infrared behavior as the one-loop scalar functions considered 
in section |l]. Finally, in appendices ^ and [b|, we collect some properties of the hypergeometric functions and the 
formulae that define the three- and the four-body phase-space in arbitrary space-time D-dimensions. 

Throughout all this paper we work with the following set of dimensionless Lorentz-invariant variables 

yij = 2(pi ■ Pj)/s, r b = mljs, 

where the are the scalar product of all the possible pairs of final state momenta normalized to s = q 2 , the 
total center of mass energy, p\ = p\ = ml is the mass of the quark/antiquark (in particular the 6-quark mass) and 
particles labeled as 3 and 4 represent the gluons, p\ = p\ = 0. Energy- momentum conservation imposes the following 
restrictions 

2/12 = 1 - 2r b - 2/13 - ?/23, 

2/12 = 1 - 2r b - 2/13 - 2/14 - 2/23 - 2/24 - 2/34- (1) 
in the case of the decays Z — > bbg and Z — > bbgg respectively. 



II. ONE-LOOP INTEGRALS CONTRIBUTING TO THREE-PARTON FINAL STATES 

We consider here the IR divergent virtual corrections to the three-parton decay q — > p\ + p-i + p% with all external 
particles on-shell, see fig |^. 

Using the Passarino-Veltman procedure we can reduce any set of vectorial or tensorial one-loop integrals to 
n-point scalar functions. Therefore only scalar integrals Jl^] have to be computed. At the NLO the relevant one-loop 
Feynman diagrams contribute only through their interference with the lowest order Born amplitudes. Therefore, we 
will be interested only in the real part of such n-point scalar one-loop integrals. Although we will not specify it, the 
following results refer, in almost all cases, just to the real part of these functions. 



A. Three-point functions 

First we present results for the two three-point functions which actually define the non-trivial IR structure of the 
whole virtual radiative corrections for the process. As we will show below, the IR poles of the four-point functions 
contributing to our process are also controlled by these simpler three-point functions. Thus, we have to consider the 
following scalar one loop integrals, see fig [lj 
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In general, this kind of functions depend on all the masses in the propagators, on the square of the difference of 
the two external momenta and on the square of the external momenta itself. Since we have fixed the mass and we 
have imposed the on-shell condition the only remaining relevant arguments of these functions are the two-momenta 
invariants j/13 and j/12, respectively. 
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FIG. 1. Infrared divergent triangle Feynman diagrams. 

For the C03 function, which has a double IR pole, the most convenient Feynman parametrization is the one with 
both Feynman parameters running from to 1 because then the two integrations decouple. After performing the 
integration over the loop momenta we have 



C03(y 13 ) 



- 1 - e /i 2e (47r)T(l + e) / x^- 2e dx 
Jo 



(4) 



x / K 1 - y)( r b - y(n + 2/13)) - irj\ 



-l-e 



dy 



We see from here that, as we mentioned, the infrared problem is at the border of the Feynman parameter space, in 
this particular case at x = and y = 1. The full calculation has to be performed for arbitrary e values. 

We will omit the ir\ term in the following calculations, however, when needed, we will use the prescription — > 
Tb — it], 2/13 — * 2/13 + if] to resolve the ambiguities. 

After integration over the x-parameter we are left with 
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(5) 



where yo = rj,/(r6 + 2/13)- The last integral gives rise to a hypergeometric function, eq. ( |A1[ ), 

C03(y 13 ) 



s -i-e^ (^Y r(i + e)r(- e ) 



r(i-e) 
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2/0 ' 
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that, after some mathematical manipulations (see appendix^) can be written in terms of a dilogarithm function, 



2^1 + 6,1; l-e;— = 
2/o 



1 

2/o 



1 



1 + 2e 2 Lz 2 — + 0(e 3 ) 
2/0 



(7) 



Above we used eq. ( |A2| ) and eq. (A3). In the region of physical interest 2/0 is smaller than one, therefore to extract 
the real parts we transform £^2(1/2/0) m terms of £^2(2/0) by using the known properties of the dilogarithm function 
Jl3[ . The ambiguities are solved by using the correct ir\ prescription. 
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The final result we obtain for the real part is, 



CQ3(y 13 ) = s- 1 -^ 



1 



1 



r(l - e) y 13 

l . n 
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-log y - Li 2 (y ) - -ry 



(8) 



Notice that for convenience we have factorized a (47r) £ /r(l — e) coefficient because this constant appears also in the 



four-body phase-space, see eq. (B4). 

The C05 function has a simple pole and already appeared in the NLO calculation of two-jet production, although 
in a different kincmatical regime. We just quote the result since its calculation is straightforward 
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r(l - e) (2/12 + 2r 6 )/3i 



- - logr fc ) logci2 - 2L(y 12 ) 



(9) 



where 



L(y 12 ) = Li 2 (ci 2 ) + 



log(l - ci 2 ) log C12 - - log C12 



and 
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B. Box integrals 

In the one-loop diagrams contributing to the NLO corrections to the three-jet decay rate of the Z-boson into massive 
quarks, see fig |2|, we encounter the following set of four-point functions, 



and 



where p 13 = p 1 + p 3 . 
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D05 D06 
FIG. 2. NLO box contributions to the three-jet decay rate of the Z-boson into massive quarks. 
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Both four-point functions are IR divergent. The integral D05 has a simple infrared pole, while DOG presents double 
infrared poles because it involves a three-gluon vertex. Before the calculation of these integrals we would like to 
discuss their IR behavior. To do that we just need to split their non on-shell propagator. 

For the D05 integral we have 



1 



(k +P13) 2 - mi 



P 2 i3 



k 2 



2k-p 13 



(Pis - m l)[( k + Pis) 2 - ml] 



(14) 



The infrared divergence is isolated in the first term of the righthand side of the previous equation that give rise to 
C05(yi2)/yi3, where C05 is the three-point function defined in eq. The other term generates two infrared finite 
one-loop integrals. 

For the D06 function things are quite more complex. The splitting of the non on-shell propagator does not give 
directly the divergent piece due to the presence of double infrared poles. First of all, we have to notice that D06 
is invariant under the interchange of particles 1 and 2, i.e., the result of the integral should be symmetric in the 
two- momenta invariants yi 3 and 1/23- The proper way for extracting the infrared piece is as follows Jl4j 



•D06(yi 3 ,y 2 3) 



Vi3 
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finite terms 
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The box integral, D05, comes from an electroweak-like Feynman diagram and it has already been calculated [|f5|] by 
using a photon mass infrared regulator. Since only simple infrared poles appear in the D05 function, we are allowed 
to identify the logarithm of the gluon (photon) mass, log A 2 , in the result of rcf. with a pole in dimensional 
regularization: log A 2 — > (47r) e /T(l — e)/e. Thus, we have 



D05(y 13 ,y 12 ) ^s^^ 
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- - log ) logci2 
e V13) 

tt 2 +log 2 c-Li 2 (l-c 2 12 ) 
2Li 2 (l - C12 c) + 2Li 2 (l - c 12 /c) 



with 



13 = VI - 4r b , c = 



1-/3 
1 + 13 



(16) 
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To calculate D06 we consider the following set of Feynman parametrizations 



1 



dz 



k 2 (k+p 3 ) 2 J [k 2 + 2k- P3 z] 2 ' 
1 

P + Pis) 2 - m l\ i( k ~ P2) 2 - m 2 } 
dy 



/ [k 2 + (2k ■ q + 2pi • P3 )y - 2k ■ p 2 ] 2 ' 
that again we combine with the help of the identity 



1 



= / dx- 



Qx(l — x) 



a 2 b 2 J [ax + b(l -x)} 4 ■ 
After these parametrizations our integral becomes linear in the z-Feynman parameter 



(18) 



(19) 
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i _ f 1 

-— — jD06(yi3, 7/23) = M 4 D / 6x(l - x)dxdy 



(2tt) d 7 [fc 2 - (M lZ + M (l - z)) + i??] 4 ' 



(20) 



where we defined 



M = s (1 - x) [r b (1 - or) - (1 - z) y(l - y) - yi3xy] , 

Mi =s(l-x) [r b (1 - 1) - (1 - x) 2/(1 - y) - y 23 x(l - y)] . (21) 

Again we will omit the irj term, and, when needed, we will use the prescription T\, — > r% — irj, yi3,23 — > 2/13,23 + *?7- 
Performing the z-integration before the integration over the loop momentum we obtain 



r^06(y 13 ,y 2 3) = s~V 2£ 



16n' 
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that results, after integration over the loop momentum, into 

£06(y 13 , y 23 ) = .s- 2 - e M 2£ (47r)T(l + e) 
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<ix (1 — x) 
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2/132/ . 



-1-e 



+(2/13 <-> 2/23) , 

where 

-4 = r b - y(l - y) . 

After integration over the x-Feynman parameter we get a hypergeometric function 

: r(i + e ) 



(23) 
(24) 
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■ dy- 
2/132/ -2/23(1 -y) 
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2/132/ n 



(25) 



Again, using the transformation properties of 2-F1 and expanding in e, we can write the above hypergeometric function 
in terms of the dilogarithm function, see appendix [a|. Moreover, as the dilogarithm function contributes at the order 
0(e 2 ) we are allowed to substitute it with its value at y = 0, £12(1) = 7r 2 /6, obtaining 



^06(y 13 ,y 2 3) = s- 2 -V e (47r) 



s r(i + e ) 
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Partly expanding in e to the required order, we can rewrite our integral before the last integration in the following 
form 



£>06(2/ 13 ,2/23) = s /" 



-2-e..2e ( 47rr &) 
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f (-2/132/)' 1 ' 26 e log(l-y(l-y)/r b ) 
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(27) 



where we have used T(l + e)T(l — e) = (1 + e 2 7r 2 /6 + 0(e 4 )). The first part of the integrand gives, up to a factor 
— (— 2/i3)~ 2e /2/i3/2/23/(— 2e) a hypergeometric function that can be reduced to the dilogarithm, 
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Using the above equation and the symmetry j/13 <-> j/23, we can rewrite D06 (1/13, 1/23) as 



^06(yi3,2/23) = s- 2 -y e #^r(l + c 2 V 
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2e 2 [ 2Li 2 



-2/23 
2/13 
2/13 + 2/23 \ 
2/23 / 
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Performing the last integral, and fully expanding in e (in this expansion the irj prescription should be used to avoid 
ambiguities) we have the following final result for the real part of the function Z?06(t/i3, 2/23) 



^06(^3,2/23) = s- 2 - e M 2e 
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III. INFRARED DIVERGENT PHASE-SPACE INTEGRALS 



We consider now the tree- level four-parton decay q — > pi + P2 + P3 + P4. When at least one gluon is radiated from 
an external quark (or gluon) line, the propagator- factors, l/(pi • Pj), can generate infrared and gluon-gluon collincar 
divergences at the border of the four-body phase-space. In this section we integrate these factors in a thin slice at 
the border of the phase-space and we show how these basic phase-space integrals are related to the infrared divergent 
scalar one-loop integrals of the previous section. 



A. Integrals containing soft gluon divergences 

The phase-space of n + 1 particles can be written as the product of the n-body phase-space times the integral over 
the energy and the solid angle of the extra particle. In arbitrary D = 4 — 2e dimensions we have 

dPS(n + 1) = — \ n _. E D - 3 dE dQ dPS{n) . (31) 
2(27r)-^ 

Suppose E3 is the energy of a soft gluon: £3 < w where w, with w very small, is an upper cut on the soft gluon 
energy. Let's consider 



/" 

J 



i^3^_l_ ; (32) 



where p$ is the momentum of the gluon and p\ is the momentum of a quark. After integration over the trivial angles 
we find (x = cos 9, with 9 being the angle between pi and p"s) 



4r(f-l)7 3 7-i (Ei-pia:) 



dE 3 E^- b / dx\- J —-r . (33) 



Here pi stands for the modulus of the threcmomcntum. In the cikonal region we can suppose E\ being almost 
independent from E 3 
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(34) 



with t/24 = 2(^2 • Pa)/ s. The integral over E 3 can be easily done and gives a simple infrared pole in e. Since the x 
dependent part is completely finite we can expand in e before the integral over this variable is performed. The final 
result we get is as follows 
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with p2 the momentum of the antiquark. In this case we can combine the denominators as follows 



(35) 
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dy-. 
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(2pi • p 3 ){2p 2 ■ p 3 ) Jo [2^3 • (pi + (P2 - Pi)y)} 



(37) 



We get therefore the same integral structure as in eq. (pq), but with the fourmomentum pa = Pi + (P2 —Pi)y, instead 
of pi . After integration over the y-parameter we get 
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2(2tt) 
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E?- 3 dE 3 dn— — 

3 (2pi • p 3 )(2p 2 ■ p 3 ) 
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Notice that this last integral has the same divergent structure as the scalar one loop C05 function defined in eq. (|9J). 
The finite contribution, the F(yi 4l y 24 ) function, is rather involved. We write it in terms of the variables yi 3 and y 23 
for later convenience 



F{yi 3 ,y 23 ) 



1 



bVT 



- T [G(z 1 ) + G(z 2 )-2G(l)} , 



(39) 



where 



G(*)=log*ilog 



h + z 



z . 1 + i z + t2 
log t 2 log ■ 



Z — t 2 



, , *i -*2 , (*i + «)(«- *a) 
+ log -^5— log - . 

2 (£i-z)(z + i 2 ) 

+ i [log 2 (i 1 +z) + log 2 (t 1 -z) 

+ log 2 (z-t 2 )-log 2 (z + t 2 )] 
1 " 



4 
iz 2 

1/12 



log(ti + z) + log(ii - 2) + log 



z - t 2 
z — tl 
z + t 2 
z + U 



Li 2 
Li 2 



z — tl 
Z + tl 
z - t 2 

z + U 



(40) 
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with 



b 2 = l- 4r b - 2(1 - 2r b )(y 13 + y 23 ) 

+(1 - r b)(yf 3 + 2/23) + 2/132/23(3 - 2r b - y 13 - y 23 ) 



where h p is the function that defines the limits of the three-body phase-space, see eq. (B3), and 

ti, 2 = (1 ± \A - a 2 )/a , 



z\ = exp 



cosh 



1 (^{l-y 13 f~Ar b /(a(l - y 13 ))) 
z 2 = exp cosh -1 U/(l - y 23 ) 2 - Ar b /{a{\ - y 23 )) 



(41) 



(42) 



B. Integrals containing gluon-gluon collinear divergences 

Let's consider the following phase-space integral 



PS{A)—^— . (43) 

2/132/34 



We work in the so-called "3-4 system", eq. (B4), where as usual p\ denotes the momentum of the quark and p 3 and 
P4 are the two gluon momenta. This integral has a singularity when either the gluon labeled as "3" is soft or the two 
gluons are collinear, that is 2/34 = 0. Then we take 2/34 < w with w very small. We can decompose the four-body 
phase-space as the product of a three-body phase-space in terms of variables 2/134 and 2/234 times the integral over 2/34 
and the two angular variables v and 9' 

P5(4)_L_ = PS(3)(y 134 ,y 23A ) 

2/132/34 

s ^ 1 r d0 > S in-^e> 



16tt 2 T(l - e) N e , J Q 
x / dv(v(l-v))-< I™ dy 3A V -^ , (44) 

Jo JO 2/13 

where S = 1/2! is the statistical factor. In the "3-4 system" the two-momenta invariant yi 3 can be written in terms 
of the integration variables as 

2/13 = ^ fj/134 - \/ 2/f 3 4 - 4 ^2/34(l - 2v)J . (45) 

The 2/13 factor is independent of the #'-angle so we can get rid of this first integral. Moreover, yi 3 contains a piece 
independent of the v- variable and another one that is odd under the interchange v «-> (1 — v). With the help of this 
symmetry we can avoid the use of the square root in eq. (|45|) by using the identity 



1/1 1 



then, we rewrite our integral as 



2\g + f(v) g-f(v)J g 2 ~f(v) 2 ' 



PS(A)^— = PS (3) (2/134, 2/234) 

2/132/34 

x if r -«))"' / 2/34^2/34 



(46) 
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22/134 



2/?34 - (2/134 - 4^2/34) (1 - 2vf 



(47) 
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After integration over 1/34 we get the first infrared pole 
PS(4) 



1 =PS(3 h S ^ ~ 1W ~ 



2/132/34 



167r 2 r(l-e) e 2y 13 4 
x / dv{v{l-v))- 1 - e 2 F 1 [l,-e-l-e--B 



(l-2^) 2 
v(l-v) 1 ' 

with B = w(rb/yf 34: ). Similarly to eq. (Eq) we can write 2F1 hi terms of a dilogarithm function 



(48) 



2 Fi[l,- e ;l-e;--] 



1 + e A Li 2 



a + b 



+ 0(e 3 ) 



Thus, we obtain 



(49) 



PS{4) 



1 



= PS(3)- 



S (47r) e -1 uT 



16tt 2 T(l - e) e 2y 134 J 
x (v(l-v) + 5(l-2i;) 2 ) e 

B(l-2i;) 2 



1 + e^i 2 



u(l - u) +.8(1 - 2d) 1 



dv(v(l — v)Y 



0(e 3 



(50) 



The above integrand is symmetric under the interchange v <-> (1 — We perform the integral only in half of the 
integration region and apply the following change of variables 



u = Av{l — v) 



(51) 



In addition, since the IR divergence occurs at u = and the term proportional to the dilogarithm function is already 
order e 2 , we can substitute the dilogarithm function by its value at u = 0, that is ££2(1) = 7r 2 /6. After this, we get 



PS(4)- 



1 



2/132/34 



PS*(3) 



167r 2 T(l - e) e 2/134 



S (l- U ) 



-1/2 



1-1 )u 

45 



(52) 



which gives again a hypergeometric function 

1 



PS (4) 



= PS(3) S (4?r)£ -1 1 2 
2/132/34 167r 2 r(l-e) e 2/134 V 2/134/ 

T(i-2e) L ' '2 ' 4£ J 



l + e' — 



(53) 



The hypergeometric function is already 1 + 0(e 2 ) but now it is not possible to write it in terms of a dilogarithm 
function. If we are interested just in the pole structure we can stop here. (Observe that neither the double pole or 
the simple one depend on the eikonal cut w). To get the complete finite part further mathematical manipulations 
must be performed with the hypergeometric function. We use |L6| 

2-Fif-e, -2e; - - 2e; 1 - — 1 = 
1 ' '2 ' AB l 

1 



{AB) 2e T(- - 2e 



2 

r(e) 



r ( - e) r(l)^ [ - 2e '2- e;1 - e;4i? ] 

iW^) 2Flh ^ ;1 + e;4S] 



(54) 
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For w small enough, w <C 2/134/^6, we can use the Gauss series [[l6| to expand the hypergeometric functions around 
B —> 0. The final result we obtain is as follows 

dPS (4)—!— = dPS(3)(y 13i ,y 2M ) " U;Ti ' ' 



2/132/34 16tH r(l - e) yi34 

_L _ ' log » „ + ( 1 _ log „) log _2_ _ 1 log * _g_ _ =_ + 0W J . (55) 



Notice we get the same infrared poles as in the one-loop three point function C03, eq. (Q), if we identify 2/134 with 2/13. 
In the limit 2/34 — > the momentum P34 = pj, + pa behaves as the momentum of a pseudo on-shell massless particle 



since P34 — ► and the function that defines the limits of the four-body phase-space, /134 in eq. (B6), reduces to the 
three-body phase-space function h p , see eq. (B3). Therefore, within this limit, this identification provides the key to 



see the cancellation of the infrared divergences. 



IV. SUMMARY 



We have presented the basic ingredients needed to compute the second order QCD corrections to the three-jet decay 
rate of the .Z-boson into massive quarks. This includes, for massive quarks, the first calculation of the double infrared 
divergent box integral, we called DOG, in dimensional rcgularization. We have presented results for both the infrared 
part and the finite contributions of some basic infrared divergent one- loop and phase-space integrals. We showed they 
share the same infrared behavior. Therefore, the cancellation of the infrared divergences is assured as stated by the 
BN and KLN theorems 01. 



We would like to acknowledge interesting discussions with S. Catani, A. Denner, A. Manorial and A. Pich. We 
are also indebted with S. Cabrera, J. Fuster and S. Marti for an enjoyable collaboration. M.B. thanks the Univ. de 
Valencia for the warm hospitality during his visit. The work of G.R. and A.S. has been supported in part by CICYT 
(Spain) under the grant AEN-96-1718 and IVEI. The work of G.R. has also been supported in part by CSIC-Fundacio 
Bancaixa. 

APPENDIX A: SOME PROPERTIES OF THE HYPERGEOMETRIC FUNCTION 

Here we present some properties of the hypergeometric function Jl6| we use in the main body of the paper. The 
integral representation of the hypergeometric function is 

2 F 1 [a,b;c;z] = / * 6 ~ 1 (1 - t) c-6-1 (± - tz)-' 1 . (Al) 



r(6)r( c -6) Jo 

From the above equation one can derive the following transformation properties, 

2 F 1 [a, b; c; z] = (1 - z) c ~ a ~ b 2 F 1 [c - a, c - b; c; z] 

= (l-z)- b 2 F 1 [b,c-a;c:-^j} . (A2) 

Using the Gauss series for the hypergeometric function m& one can easily see that 



2 



Fi[ae,/?e; 1 + je; z] = 1 + e 2 af3Li 2 (z) + (9(e J ) , (A3) 



where, Li 2 {z) is the Spence or dilogarithm function 



Ll2[z) = I dt l ^-A . ( A4 ) 
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APPENDIX B: PHASE-SPACE IN D = 4 - 2e DIMENSIONS 



The phase-space for n-particles in the final state in arbitrary D space-time dimensions [0 (D = 4 — 2e) has the 
following general form 

dPS{n) = (2n) D 11 d °~ lp 



= 1,71 

d D P, 




- (2 " D n (2lt) o-. 

x 3(Pi ~ rrii)<d(E l )6 D I </ > I . 1 131) 




Let's consider the decay into three particles, q — > p\ + pi + P3, where particles 1 and 2 share the same mass, 
p\ = p\ = m^, and particle 3 is massless, p\ = 0. In terms of the two-momenta invariant variables j/13 = 2(p\ -pz)j s 
and j/23 = 2(p 2 ■ 2>3)/», with s = q 2 , we get 

P5 ^ 3 ^ = 16(2tt) 3 T(2 - 2e) (iff) 

x /" 9(h p )h; e dy 13 dy 23 , (B2) 

where the function /i p which gives the phase-space boundary has the form 

h P = 2/132/23(1 - 2/13 - 2/23) - H,(yi3 + to) 2 , (B3) 

with ?*b = rn^/s. 

For the case of the decay into two massive and two massless particles, q — ► pi + P2 + P3 + £>4 with = p| = TO b 
and P3 = p\ = 0, it is convenient to write the four-body phase-space as a quasi three-body decay 

Q -> P34 + Pi + P2 
^ P3 + J?4 • 

In the cm. frame of particles 3 and 4 the four-momenta can be written as 

pi = (£1,..., 0,pi) , 
P2 = (Eh, ■ ■ . ,p 2 sin^, p 2 cos/3) , 
P3 = Es(l, . . . , sin 6* cos 6*', cos 0) , 
Pa = E±(l, . . . , — sin 8 cos 6', — cos 0) , 

where the dots in p$ and P4 indicate Z? — 3 unspecified, equal and opposite angles (in D dimensions) and D — 3 zeros 
in pi and p 2 ■ We will refer to this as the "3-4 system" M . 

Energies and threemomenta can be written in terms of the following invariants 

2fes ' Pi) 

2/34 = > 

s 

2(P1"P34) 2(P2-P3i) 

2/134 = , 2/234 — , 

S S 

where p^ = P3 + pa- We obtain 

P 2/l34\/s v 7 ^ f~2 'a 

e i = o , > Pi = .— V 2/134 - 4^2/34 , 

2^2/34 2^yil V 

„ 2/234\/s v 7 * /~1 , 

^2 = ; , P2 = 7TT= V ^234 ~ 4r &2/34 , 

2Vy34 2^/^37 v 



£/3 — £/4 — 
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Setting v = (1 — cos#)/2, the £)-dimensional phase-space in this system is 

PS(4)= 16(2*)' r(2-2e) (^) 2 7 dWy234 

„l-e ^ (4tt) £ 



XS £ 16tt 2 r(l-e) ' dy34d( - h ^ h 3i y 3i 



sin - * 9' , (B4) 



where S = 1/2! is the statistical factor, Ng> is a normalization factor 



V, ; ,= / ^ Sin - 2 ^' = 2 2 %^L^1, (B5) 



and the function 



hu = 2/1342/234(1 - 2/134 - J/234) ~ n>(2/l34 + V23i) 2 

+ (-1 + 4r b + 2(1 - 2r b )(yi 34 + 2/234) 

" 2/134 " 32/1342/234 ~ 2/234) 2/34 

+ 2(1 - 2r b - yi 34 - 2/234)2/34 ~ 2/34 » ( B6 ) 



defines the limits of the phase-space. Observe that /134 reduces to eq. (|B3j) in the case 2/34 — > 0. Furthermore, in this 
limit, p 3 4 = p 3 +P4 behaves as the momentum of a pseudo on-shell massless particle since p| 4 — > 0. 
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